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We study the segmental dynamics of poly (ethylene oxide) (PEO) from microscopic simulations 
in the neat polymer and a polymer electrolyte (PEO/LiBF4) by analyzing the normal modes. We 
verify the applicability of the Rouse theory, specifically for the polymer electrolyte where dynamic 
heterogeneities, arising from cation-polymer interactions, alter the mobility non-uniformly along the 
chains. We find that the Rouse modes for both the systems are orthogonal despite the presence of 
non-exponential relaxation of the modes and violation of the Gaussian self-similarity of the chains. 
The slowdown of the segmental dynamics in the polymer electrolyte is rationalised by an order of 
magnitude increase in the friction coefficient for those monomers which are bound by cations. In 
general, for the electrolyte the Rouse predictions for the dynamics of segments (both free and/or 
bound) agree well except for very short times. 
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Introduction 



Polymer electrolytes^^ are obtained by dissolution of 
a salt in a polymer matrix. The multitude of interactions 
occurring in these systems gives rise to complex dynam- 
ical processes of the ions and the polymer segments over 
a broad range of length and time scales. On the one 
hand, comprehending the relevance of the mechanisms of 
ion dynamics will allow the prediction of ionic mobilities, 
and, on the other hand, the underlying polymer dynam- 
ics attains significance due to the coupled nature of the 
motion of ions and polymer segments. 

Computer simulations have provided precious insights 
on the structural aspects as well as the dynamical pro- 
cesses in polymer electrolytes. S 01 IS El Hi The cations of 
the salt are solvated by atoms or groups located on the 
polymer host by means of complexation. This results 
in a local increase of the friction coefficient of the com- 
plexed monomers whereas the uncomplexed monomers 
are relatively free. The difference in the dynamics of 
the complexed and the uncomplexed monomers lead to 
dynamic heterogeneities.^" Molecular dynamics (MD) 
simulations^ of lithium salt in poly(ethylene oxide) 
(PEO) have shown that the dynamics of oxygen atoms 
which are bound to a lithium cation (Li + ) are slower 
than those which arc free.M The mean square displace- 
ments (MSD) of oxygen atoms exhibit a time exponent 
of t a with a w 0.6 for intermediate times less than the 
relaxation time of the end to end vector of the polymer. 

The widely studied Rouse theory,^ El known to ap- 
proximately describe the dynamics of linear, unentangled 
polymers under melt conditions, capture the essence of 
the dynamical behavior of the segments. The study of 
polymer dynamics can be broken down into several dy- 
namical modes with characteristic time and length scales. We have performed MD simulations of a system of i) 

neat PEO and ii) a model polymer electrolyte composed 

of PEO and a salt, LiBF4. The concentration of the 

polymer electrolyte is 20:1 in terms of the EO:Li + ra- 

* Electronic address: |arijitmaitra@uni- muenster.de, andheuerQuni^mu erfealiar.del'he GromaCsE3 simulation package was used for 



Within the framework of the Rouse theory, each mode 
corresponds to an eigenvector (also called normal coordi- 
nate). The eigenvectors are orthogonal to each other, i.e., 
the different modes are statistically uncorrelated. Knowl- 
edge of modes can be employed to calculate the segmental 
dynamics. The Rouse theory predicts for the segmental 
MSD a dependence of t - 5 for times t -C £r where tn 
is the Rouse time or the longest relaxation time of the 
polymer. 

Very recently, it has been shown for PEO/LiBF4 that 
it is possible to predict the cationic diffusivity as a func- 
tion of the microscopic time scales characterising the 
cation transport mechanisms by considering the coupling 
of the cationic motion to the segmental dynamics of the 
polymer. EB An important delineation was the descrip- 
tion of the ether oxygen (EO) dynamics according to the 
Rouse theory, which was verified by comparison against 
the MSD of EO atoms, and particularly to account for 
the retardation of those monomers which are bound to 
the Li + ions. 

In this paper we elucidate the polymer dynamics 
from microscopic simulations of amorphous PEO and 
PEO /LiBF4 and check the extent to which the dynamics 
can indeed be characterized by the Rouse model. Fur- 
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thermore, we present Brownian dynamics 1 simulation 
of a bead-spring model with heterogeneous friction and 
some stiffness to closely resemble our chemically realistic 
polymer electrolyte. In this way we analyze to which de- 
gree the dynamic properties of the microscopically com- 
plex PEO /LiBF4 system can be understood from an ap- 
propriately chosen bead-spring model. 



Models and Methods 
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generating the MD trajectories. All the atoms — viz. 
carbon, oxygen, hydrogen, lithium, boron and fluorine 
— were explicitly taken into consideration for setting up 
the respective systems. The nonbonded interactions be- 
tween all pairs of atoms including those of the cations 
and anions were modeled according to the Buckingham 
potential plus the electrostatic interactions. The two- 
body effective polarizable force fields employed have been 
adopted from ref . ^ ^ which were derived from quantum 
chemistry based techniques. The simulation methodol- 
ogy and details of the force fields used in this work closely 
follows that of ref. SC1 

Coulomb interactions between the partial charges (ions 
as well as polymers) were calculated via the particle-mesh 
Ewald method with a distance cut-off of fO A. Overall 
charge neutrality of the systems were also maintained. 

The simulation box contained 16 chains of PEO, 
with the sequence H-(CH2-0-CH2)at-H, each comprising 
N = 48 number of monomers. All bond lengths were 
constrained to their equilibrium values by the L1NCS 
algorithm^. However, the bend and the dihedral an- 
gles of the polymer chains were kept flexible. 

The simulations were propagated at a constant tem- 
perature of 450 K in the canonical ensemble and under 
the application of periodic boundary condition. Temper- 
ature coupling was accomplished using the Nose-Hoover 
thermostat.^ ^ The average pressure during the runs 
were of the order of f MPa. Equilibration times of the 
order of tr had been considered for each of the systems 
prior to the analysis of the trajectories. 

Apart from the realistic systems described above, we 
have also simulated a toy model of a single chain con- 
sisting of only beads and springs, characterized by the 
standard harmonic potential. We attempt to capture the 
essential dynamics of the monomers of our realistic poly- 
mers through the toy model. To incorporate correlations 
between the bond vectors we impose a bending potential 
such that the total Hamiltonian is written as 



h N ~ 2 

N-2 
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where is the position vector of bead i; k and kg are the 
force constant of the harmonic springs and the bending 
constant between the adjacent bond vectors, respectively. 
The forces arising from the bending potential essentially 
model a freely rotating chain (also see ^). A value of 
k = 3ksT is chosen, where ks — 1 is the Boltzmann 
constant in reduced units and T = 1 is the temperature, 
also in reduced units, used for coupling between the vis- 
cous drag and the Gaussian white noise, that entails the 
mean square bond length I 2 = 1.0. We fix kg = 2ksT to 
maintain consistency with the stiffness of the PEO chains 
in the atomistic PEO/LiBF4 system, yielding a similar 
value of the characteristic ratio Coo for both the models; 
see Table f . 



We construct three different toy systems with system 
size N = 48: (i) BD1: the chain is homogeneous, i.e. all 
the beads have the same friction coefficient of Ci = lj 
and akin to the characteristics of neat PEO. (ii) BDRA: 
12 beads are chosen randomly with the friction coeffi- 
cient C2 = 8 to reflect heterogeneity (or defects), (iii) 
BDCO: two sets of 6 contiguous defects are chosen such 
that the distance between any two particles belonging to 
different sets is at least 6 particles apart. This realiza- 
tion qualitatively resembles contiguous monomers (« 5-7 
in number^) of PEO bound by cations as in the 20:1 
PEO /LiBF4 system. The choice of the friction coeficient 
of the defect £ = 8 has been rationalised later in this 
article. 

The Brownian dynamics method is employed to evolve 
the toy systems with a time step of 0.001. For data anal- 
ysis, all the computed quantities from the toy systems 
are averaged over 20 independent runs. For BDCO the 
distance between the two sets of contiguous defects is 
chosen randomly for each run. 



Chain Static Characteristics 

The global chain properties are listed in Table 1. The 
correlation between adjacent bond vectors measured as 
(cos 9) where 7r — 9 is the angle included between the two 
adjacent bond vectors, r^— r,_i and r.; + i — r,, is consistent 
between the atomistic systems and the toy models. For 
PEO, 7r — 9 is the angle between the link vectors formed 
by connecting the ether oxygen atoms consecutively. The 
stiffness of the PEO chains as measured by (cos 9) is 
slightly higher in the neat system when compared to 
PEO/LiBF.4. This implies the chains are expanded more 
for the neat PEO than those in the PEO/LiBF4 system. 

Alternatively, this is also apparent from the values of 
the characteristic ratio, C oa =R 2 /l 2 (N - 1),^ where R 2 
is the mean square end to end distance of the polymer 
chain, I 2 is the mean square bond length and TV — 1 is 
the number of bonds. Note, that in ref. ^ Coo ~ 5.1, 
for the case of neat PEO, because I 2 was defined as the 
mean square length of a skeletal bond. We consider I 2 to 
be the mean square O — O distance between the adjacent 
monomers. The characteristic ratio is not sufficient to 
describe the chain stiffness because it takes into account 
only the large scale chain conformations. Even though 
Coo is higher for neat PEO when compared to the chains 
in PEO/LiBF 4 , the chains in the former display faster re- 
laxation (tr ki 10 ns) than those in the electrolyte (tr ss 
19 ns). This underscores the importance of local stiff- 
ness effects arising from dynamic heterogeneities which 
are difficult to incorporate in analytical models. Struc- 
turally, as evident from the O — O radial distribution 
function (not shown) and also from the mean square dis- 
tance between adjacent oxygen atoms, l 2 , there will be 
slight differences between the PEO chains in the two sys- 
tems in terms of the distribution of the conformational 
states of the -C-O- and -C-C- dihedrals. 13 ^ For 
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System 


i 




^COS f7 / 


PEO 


10.0 [A 2 ] 


3.3 


-0.63 


PEO/LiBF 4 


9.2 [A 2 ] 


2.7 


-0.55 


BDl,BDRA,BDCO 


1.0 


2.6 


-0.53 



Table I: Chain static properties of the different model sys- 
tems from simulations. The units for the microscopic simu- 
lations (PEO and PEO/LiBF4) are provided in box brackets 
and data for the toy models (BD1, BDRA, BDCO) are ex- 
pressed in reduced units. (I 2 : mean square bond length be- 
tween monomers/segments; R 2 : mean square end-end vector 



of the chains; characteristic ratio C a 
mean of the cosine of bond angles) . 



R 2 /l 2 (N-l); {cos9): 



the toy models all the static properties, which are inde- 
pendent of friction coefficients of the beads, are similar 
as is expected. 



Analysis of Rouse Modes 

Statics 

The static and dynamic properties of the polymer 
chains can be studied by means of the normal mode anal- 
ysis. The normal modes (also called Rouse modes) for a 
discrete polymer chain^ is expressed as 



1 



JV 



n=l 



(n — 0.5)pir 
N 



(2) 



where p = 0, . . . , N — 1 and r„ denotes the center of mass 
of monomer n. The static mode amplitudes, (X 2 ) of the 
different models are shown in Figure 1(a). A dependence 
of sin~ 2 (piT/2N) w p~ 2 would conform to the Rouse pre- 
diction which essentially considers the distribution of the 
normal coordinates to follow the Gaussian distribution. 
However, owing to the chain stiffness systematic devia- 
tion is observed for all the models studied (see also ^). 
All the Brownian systems, BD1, BDRA and BDCO show 
the same static p-behavior (not shown). The homoge- 
neous Brownian chain (BD1) shows a dependence of p~ 22 
only for the first few modes and then changes to a steeper 
exponent reported earlier. The Brownian chain systems 
behave effectively as a freely rotating chain (FRC). The 



analytical result 
given by 



24] 



for the normal mode amplitudes is 



with 



R 2 



8N(N - 1) 



(l + (cosfl)) ; 
4(cos6>) 



2 p^y 1 

(3) 
(4) 



A fit of this equation with (cos 9) = —0.45 agrees well 
with the simulation data (see Figure 1(a)). This angle is 
close to the actually observed angle of (cos#) = —0.53. 

The neat PEO system shows a p~ 2 dependence for 
small modes up to p — 5. Beyond this the relation- 
ship smoothly curves to an approximate dependence of 
p~ 4 . Similar violations have been reported before for 
1,4-polybutadiene^" and n C ' c.:: 1 1^:;^ - However, the 
p-dependence for the high frequency modes is much 
stronger in case of PEO. The self-similarity of the Gaus- 
sian coil over a broad length scale is no longer fulfilled. 
The cause can be attributed to the preference of se- 
quences of conformational states that are more coiled 
than others.^ 

The polymer electrolyte system also exhibits a p~ 2 de- 
pendence for only the first few modes and crossing over 
sharply beyond p « 20 to a p~ 5 dependence. Specially, 
for high modes the presence of cations in the proximity 
of some monomers leads to a localised shrinking of the 
affected subparts of the chains. Precisely, the cations en- 
force a helical like arrangement of the ether oxygen atoms 
around themselves and thereby further reducing the ex- 
tent of spatial occupation of a contiguous group of bound 
monomers. This is consistent with the reduction in R 2 
and also supported by the absence of such observations 
for the highest modes in the BD1 system where there is 
no torsional potential and non-bonded interactions. 



Dynamics 

Not surprisingly, the models BDRA and BDCO show 
biexponential decay of the autocorrelation functions of 
the modes because of the presence of two different kinds 
of particles in a chain. This effect becomes very pro- 
nounced with increasing mode index. A similar effect 
is expected for the PEO/L1BF4 system, albeit diluted, 
due to frequent making and breaking of linkages with the 
ether oxygen atoms and intrachain cation motion. El 03 
non-exponential fit, like the Kohlrausch- William- Watts 
(KWW) fit, nevertheless, yields fair agreement for the 
microscopic models and BD1. 

The Rouse theory predicts an exponential decay of the 
autocorelation function of the p-th mode: 



(X p (<)X p (0)) = (X2)cxp(-i/r p ). 



(5) 



where r p is the time constant for the p-th mode. In gen- 
eral, deviation from the exponential behavior has been 



observed before. 



mm 



We have extracted the parameters 



of the KWW form, exp{ — (t/T p )^}, where the stretching 
parameter, (3, characterises the non-exponentiality of the 
decay process. Figure 1(b) shows the relationship be- 
tween Tp/ri and the mode index p for the different sys- 
tems. The Rouse theory prediction of t p oc p 2 is observed 
only for the first few modes for the systems, reminiscent 
to the observation of the p-dependence of mode ampli- 
tudes. 
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Figure 1: Analysis of Rouse modes. (Note that the data sets 
from the different models are offset for clarity in viewing) (a) 
Rouse mode amplitudes (X?) scaled by the amplitude (X?) 
of mode p=l as a function of the mode index p. A fit using 
the Equation @ for a freely rotating chain with {cos 8) = 
—0.45 is shown by the dotted line, (b) Scaled relaxation time 
of pth mode, t p by t\ in dependence of mode index p. (c) 
Stretching parameter (3 from the Kohlrausch- William- Watts 
(KWW) fitting, exp(— t/r p )' 3 , to the mode autocorrelation 
function (X p (t)X p (0)} as a function of p. 



Figure 1(c) shows the relationship between (3 and mode 
index p. For the Brownian chain, BD1, the (3 values 
fluctuate around 0.9. Thus, the correlation functions 
for all modes are nearly exponential. The PEO and 
PEO /L1BF4 systems show consistently decreasing values 
(3 with increasing mode index. The non-exponentiality 
of the relaxation processes progressively increases with 
decreasing length scales due to the existence of several 
structural realizations of a sequence of monomers in the 
presence and/or absence of a nearby cation. In gen- 
eral, the non-exponentiality of the relaxatory modes is 
stronger for the electrolyte than for the neat polymer 
due to the existence of cations in the vicinity of some 
monomers in the former. 

The Rouse theory provides a linear relationship be- 
tween (Xp) / (X^) and t p /t\. This is satisfied for the BD1 
chain as shown in Figure 2 for all the modes. The agree- 
ment in case of PEO and PEO /LiBF4 is only up to mode 
p = 12 and p = 7, respectively, and beyond which devi- 
ations are observed. The data suggest that these devia- 
tions may be related to the observed non-exponentiality 
in {X p (t)Xp(0)) which are beyond the scope of the Rouse 




X Ac, 
p 1 

Figure 2: Scaled mode amplitude {X.p)/{X.f} versus scaled 
relaxation time r p /ri of the different modes obtained for 
PEO/LiBF 4 (circle), PEO (square) and BD1 models. The 
Rouse expectation (X^) /(X?)=t /ti is shown by the gray 
lin< 
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Figure 3: Mean square displacement (MSD) of ether oxy- 
gen against time for PEO/LIBF4 (circle) and PEO (square). 
MSD calculated from Equation §6§ using the autocorrelation 
data extracted from simulations (dash). Rouse prediction of 
segmental dynamics (gray solid line). 

theory. 

Mean Square Displacement 

The mean square displacement (MSD), g(t) of a poly- 
mer segment relative to the center of mass (c.o.m) of the 
polymer chain can be expressed in terms of its normal 
coordinates: 

N-l 

9 (t)=4£[(x2)-(X p (f)X p (Q))] (6) 
P =i 

The only assumption entering Equation ((6|) is the or- 
thogonality of the modes. Figure 3 shows the MSD of 
ether oxygen atoms of the PEO and PEO /LiBF4 system. 
During the intermediate times, the segments experience 
sub-diffusional dynamics, t a , with a ~ 0.6 arising from 
the topological constraints. Using Equation © we have 
evaluated the MSD from the autocorrelation function of 
the modes. The excellent agreement between the two, for 
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Figure 4: Mean square displacement (MSD) relative to chain 
center of mass: (i) all segments (solid) and (ii) defects 
(dash) computed for the toy models BDRA (gray) and BDCO 
(black), respectively. MSD of segments for the system BDl 
are shown by dotted line. Also marked are the time scalings 
at short and intermediate times that is required for agreement 
between the relevant parts of g(t) (see text for details). In- 
set: MSD of ether oxygen (EO) of neat PEO (dots); all EO 
atoms of PEO/LiBF4 (solid) and EO atoms bound by Li + 
during time t (dash); X and Y axes are scaled for the sake of 
comparison against the toy models. 



both the systems, is indeed consistent with the orthogo- 
nality of the modes. We have also verified independently 
from the cross correlations (not shown) between the dif- 
ferent modes that the modes are orthogonal to each other 
in agreement with similar observations before.^ The 
MSD from the ideal Rouse theory obtained by inserting 
Equation ([5]) into Equation ([6]) and using t p = Ti/p 2 is 
also compared in Figure 3. For the PEO system the the- 
ory and MD simulation results agree excellently beyond 
100 ps. For shorter times, one sees only slight deviations. 
For the polymer electrolyte, beyond 20 ps the agreement 
with the Rouse prediction is very good. For times shorter 
than 20 ps the microscopic details, as captured by the 
non-exponentiality and the p-scaling of (X^), come into 
play and digressions from the ideal Rouse prediction be- 
come visible. 

An ether oxygen atom is considered to be bound to a 
Li + if their separation is less than 3.0 A. We attempt to 
elucidate the difference in dynamics observed for bound 
and average (i.e. both free and bound) oxygen atoms 
(Figure 4) through the simplified Brownian chain models 
BDRA and BDCO. 

The MSD of the average segments and the defects 
with respect to the center of mass of the polymer chains 
are shown in Figure 4. For the BDRA as well as the 
BDCO, at short times, when the particles are yet to feel 
the topological constraints, the difference in the dynam- 
ics can be expressed by the statistical factor: g(all seg- 
ments) /g(defect) = C2(0.75/Ci +0.25/C 2 ) ~ 6.3. This is 
in agreement with the data; see Figure 4. 

At intermediate times, the average segment of BDRA 



is seen to lag behind that of the BDCO whereas the defect 
dynamics is slightly faster in case of BDRA compared to 
BDCO. Furthermore, the dynamics of the average seg- 
ment and the defects merge with each other, correspond- 
ing to the value of g(t) « R 2 /3, earlier for the BDRA 
than their counterparts for the BDCO. The reason is 
that for BDRA the dynamics of the small-p Rouse modes 
quickly gets averaged over the fast and slow monomers. 
Slow and fast segments exchange their dynamics leading 
to retardation of the fast ones and acceleration of the 
slow beads until all segments share the same dynamics. 
This process occurs at a faster rate for BDRA as the 
possibility of each defect, on an average, to interact with 
non-defects is considerably higher than if the defects are 
spatially correlated. 

In contrast, for BDCO there exist extended regions 
of slow monomers which even at intermediate times can 
continue to move slower than the rest of the chain. Be- 
cause of slower exchange of dynamics between the slow 
and fast segments, the process of averaging the segmental 
motion takes a longer time compared to the BDRA. 

Additionally, we plot the MSD of the beads from the 
homogeneous chain (BDl). As expected, these show 
faster motion compared to the heterogeneous chains. At 
intermediate times, a time scaling of the MSD of the BDl 
segments by a factor corresponding to the average friction 
coefficient of the heterogeneous chain 0.75£i + 0.25(2 ~ 
2.75 (2.8 in Figure 4) is found to describe the defect dy- 
namics of the BDRA chain fairly well. 

The dynamics of the EO atoms in the neat PEO is 
displayed in the inset of Figure 4. Also shown for the 
PEO/L1BF4 system are the MSD of the average EO 
atoms and a subset of EO atoms that are bound to a 
Li + ion during time t. Note that the x and the y axes of 
all the PEO-based curves are scaled to match the units of 
the BD systems to facilitate a comparison. The oxygen 
dynamics in the neat PEO is significantly faster when 
compared to the average oxygen and the bound oxygen 
atoms of the polymer electrolyte. It can be observed 
that between t = 1 and t — 100 both the curves from 
the PEO /LIBF4 system basically agree except for a time 
scaling with a factor of 2 (see inset of Figure 4). This 
time scaling signifies the relative immobilization of the 
bound oxygen atoms in the electrolyte.^ A similar fac- 
tor of » 2 is observed in the BDCO toy system between 
the dynamics of the average segment and the defects in 
the time range of t = 50 and t = 200. Thus, a represen- 
tation of the dynamics of the complexed monomers (or 
defects) necessitates an enhancement of the friction coef- 
ficients by one order of magnitude (£ 2 /(i=8) to generate 
a right-shift of the dynamics of the average segments by 
a factor of rj 2. This is also the reason behind choosing 
a value of (2 = 8 as the friction coefficient for the defects 
in our toy systems. 

At intermediate times a time scaling by a factor of 
3 relates the EO dynamics in the neat PEO to the 
PEO/L1BF4 system. In contrast, a comparison of the av- 
erage segmental dynamics of the BDl and BDCO yields 
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a factor of 1.4. This indicates that a rigorous comparison 
of EO dynamics in the neat PEO and the PEO/LiBF 4 
is somewhat hampered owing to the differences in struc- 
ture as well as the mean square end to end distance (of 
« 20%) between the two systems. Nevertheless, the dy- 
namics observed in BDCO is in qualitative agreement to 
the dynamics exhibited in the polymer electrolyte sys- 
tem. 



Conclusions 

The polymer dynamics from all-atom simulations of 
PEO and PEO/LiBF4 are expressed in terms of its 
normal coordinates or modes. The orthogonality of the 
modes holds and is supported by the accurate repro- 
duction of the segmental dynamics of the monomers by 
taking into account only the autocorrelation functions 
of the modes and ignoring the presence of cross correla- 
tions. We have compared the findings to the predictions 
from the Rouse theory. The theoretical segmental 
dynamics considering only the exponential relaxation 
of the mode-autocorrelation function agrees with the 
monomer dynamics in case of PEO/LiBF4 except for 
the first 20 ps where the displacements are shorter 
than the statistical segment length. The agreement 
between theory and simulation holds for PEO only 
beyond 100 ps. This is surprising because one would 
have expected stronger deviations in case of polymer 
electrolytes than for neat PEO due to additional gener- 
ation of heterogeneities in the former by the complexed 
cations. For both PEO and the PEO/LiBF 4 the mode 
amplitudes show a dependence of p~ 2 only for the first 



few modes before crossing over to a steeper dependence 
for high modes where local chemical details start to play, 
overriding the Gaussian self-similarity seen over larger 
length scales. Furthermore, both the systems suffer 
from rising non-exponcntiality with increasing mode 
numbers due to diversity in local environments and 
arrangement of conformational states. The violations of 
both (Xp) and r p from the p~ 2 dependence is the cause 
of the steeper time exponents at intermediate times of 
the segmental dynamics compared to that predicted 
from Rouse theory. Finally, the effect of the presence of 
cations in the vicinity of segments augments the friction 
coefficients of the affected subparts of chain by an order 
of magnitude. Specifically, the dynamical effects of the 
heterogeneities when they are confined along the chain 
in the form of contiguous units is more drastic than 
when they are distributed over the chain. Nevertheless, 
the Rouse model captures the average behavior of the 
segmental dynamics in PEO/L1BF4 quite well. The 
dynamics of the bound segments can be obtained via 
simple time-scaling in the time regime which is relevant 
for the understanding of the ion dynamics in polymer 
electrolytes.^ 
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